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Complex poles of the unitarized tt-tt scattering amplitude in nuclear matter are studied. Partial 
restoration of chiral symmetry is modeled by the decrease of in-medium pion decay constant /*. 
For large chiral restoration {fn/ fn <C 1), 2nd sheet poles in the scalar (cr) and the vector (p) mesons 
are both dictated by the Lambert W function and show universal softening as /* decreases. In- 
medium TT-TT cross section receives substantial contribution from the soft mode and exhibits a large 
enhancement in low-energy region. Fate of this universality for small chiral restoration {f^/ f-w ~ 1) 
is also discussed. 
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An analogy between in-medium hadrons in QCD 
and collective modes in condensed matters suggests 
that hadronic spectral functions and complex hadronic- 
poles have important information on chiral structure 
of hot/dense matter [0. In particular, the light vec- 
tor mesons (p, ^,0) (2|7 and the light scalar meson (cr) 
1^ have been proposed to be possible probes of partial 
restoration of chiral symmetry. Properties of in-medium 
vector mesons may be extracted through dilepton emis- 
sion from hot and/or dense matter. In fact, relativistic 
heavy ion experiments in SPS at CERN indicate some 
spectral broadening and/or shift of the p jH [H). Dilep- 
tons observed in proton-nucleus reactions at KEK also 
suggest a low-mass spectral enhancement in the vector 
channel 

The fj meson, on the other hand, has been thought to 
be a direct but difficult probe of chiral restoration be- 
cause of its obscure nature in the vacuum. Nevertheless, 
there arise growing interests in recent years to the a in 
nuclear matter inspired by new experimental data of the 
TT-TT invariant mass distribution in (tt, tttt) [R ^ and 
(7,7r7r) reactions with nuclear targets |l^. They 
suggest some spectral change of the two-pion final states 
in the I—J—0 channel (/ and J being the total isospin 
and total angular momentum, respectively). Theoreti- 
cally, the a pole located deep in the complex energy plane 
may move toward the real-axis Q and behave as an im- 
portant precursory mode if the partial chiral restoration 
occurs [ll H, 0. 

The main aim of this paper is to investigate a com- 
mon mechanism dictating the behavior of both a and p 
when the partial chiral restoration takes place in nuclear 
matter. For this purpose, we study complex poles of the 
TT-TT scattering amplitude in nuclear matter calculated in 
unitarized chiral models. To make reliable unitarization 
and analytic continuation of the scattering amplitude to 
complex energies, we adopt the N/D method ^ |l^ 
which respects the analyticity and approximately satis- 
fies the crossing symmetry. As for the chiral models to be 
unitarized, we adopt four complementary models shown 
below. They are useful to check the model dependence 
of the final results. 



Model A: The "p model" |l9[ where tt and bare-p are 
the basic degrees of freedom. The physical a is generated 
dynamically in this model. 

Model B: The "cr model" |o) where tt and bare-cr are 
the basic degrees of freedom. The physical p is generated 
dynamically in this model. 

Model C: The "degenerate p-cr model" where tt, bare- 
p and bare-cr are all basic degrees of freedom. 
Model D: The leading order chiral lagrangian £2, by 
which the cr-pole is generated dynamically while the 
physical p is not Q. 

To incorporate the effect of chiral symmetry restora- 
tion in these models, we replace the pion decay constant 
by an in-medium decay constant /* which is expected 
to decrease in nuclear matter In the cr-models, this 
replacement is approximately justified in the mean-field 
level where the effect of the Fermi-sea is absorbed in the 
redefinition of /^r Q. In general, there arises two de- 
cay constants in the medium: /* ^ which is related to the 
temporal component of the axial current, and /* j, which 
is related to the spatial component of the axial current. 
In the qualitative study given below, we will not consider 
this cornplication. For more discussions on /*, see recent 
papers P2| . 

As for the parameters such as the masses of bare-p 
and bare-cr and their couplings to tt, we keep them den- 
sity independent, partly because it is a simplest possible 
choice and partly because their density dependence is 
not known. Nevertheless, we shall show that, in all the 
above models, there arise simultaneous softening of cr and 
p driven by the decrease of /*. This is the first strong 
indication that both the scalar and vector mesons serve 
as chiral soft modes. To extract essential physics of the 
soft modes without complication due to the finite pion 
mass, we will work in the chiral limit below. 

The TT-TT scattering amplitude and associate phase shift 
in the vacuum MeV) in Model C has been stud- 

ied in by using the N/D method. We briefiy re- 
capitulate its essential parts which are necessary for 
the analyses in this paper. The invariant amplitude 
for the elastic process tt°- -\- tt^ — > tt'^ -f tt'' is written 
as Mabcdis,t) = A{s,t)Sab6cd+ (crossing terms), where 
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a,b,-- ■ denote isospin indices and s,t are the Mandel- 
stam variables. The tree-level amplitude yl'''°°(s,t) is 
parameterized by the contact tt-tt interaction from £2 
together with the exchange of bare-cr in the s-channel 
and the bare-p exchange in t, u-channcls [ p^ . Then the 
partial wave amplitude , has a general form 



(1) 



Here m„(^p-^ is the mass of the bare-fT(-p). The first term in 
the right hand side of eq.(0) is a model-independent part 
where b^j is the low energy constant solely determined 
by chiral symmetry; 
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The second term in is a model- dependent part where 
the coefficient b'^^^ is proportional to with being a 
strength of the coupling of a with two pions. fa{s/m^) 
behaves as 0{s^) for s ^ 0. We take fhc— 00 {Model A), 
rhp=oo{Model B), and rri„,p=oo in Model D. For Model C, 



we take the same constraints as ref. 
and fhcr = fa p. 



namely, ga=gp 



In the N/ Z?-method, the full amplitude is written 
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where Njj (Dij) has a left- (right-) hand cut in the com- 
plex s-plane. The elastic unitarity implies lmDij{s > 
0) = —Njj. For the numerator function, we take 
Njj{s)=afj^''{s) together with Djj{0) = 1, which are 
consistent with the boundary condition 0^-^(5 ^ 0) ^ 
bjjS obtained from chiral symmetry. Then the disper- 
sion relation for Djj{s) reads 



Djj{s) - - 

TT 



ds'- 



(subtraction). (4) 



Since a*''™ 



(s — > 00) cx s, two subtractions are necessary 
to make Djj[s) finite. Namely two unknown parameters 
appear. One of them can be fixed by Djj(0) = 1 men- 
tioned above. Another one, which is written as d'^^^fi) 
below, should be determined empirically at some energy 
scale fi. Thus the final form of the denominator function 
reads 



Djj{s) = 




h'^,M'K) -(5) 



Here da{x) is obtained by the dispersion integral of fa{x) 
in Note that Djj{s) is ^ independent as a whole. 

For our purpose, we need (i) to find complex poles of 
ajj{s) (or equivalently the complex zeros of Djj{s)) in 
the 2nd Riemann sheet of the s-plane, and (ii) to find 
trajectories of those poles as a function of /*. In eq.(||). 



TABLE I: The bare mass and d'n extracted from the 
global fit of the phase shift in the I—J=l channel up to 1 
GeV m. For Modd^D, we use Ll{fi) = 0.4 x 10"^ and 



L^(/i) = 1.4 X 10" 



to extract d' at jj, = 0.77 GeV. doc 



obtained from d'n is also shown. 



Model 


ma (GeV) 


d'ii(GeV-2) 


d[,o(GeV-2) 


A 


mp=0.774 


-0.385 


0.171 


B 


ma=0.838 


1.66 


-3.04 


C 


mp,a=0.778 


-0.342 


-0.0855 


D 




0.277 


2.19 



log(— z) and da{z) have multi- valued structure. da{z) 
contains both ln(— z) and the Spence function Sp(z) — 
— dy ln(l — y)/y [|l6[ . The former (latter) has a branch 
cut along the real axis for Rez > (Rez > 1). Then the 
analytic continuation of the Spence function to the 2nd 
sheet reads, 



Sp"(z) = Sp^(z) + 2Tii ■ Ln(z), 



(6) 



where z is located in the lower half plane and Ln(z) = 
In |z| -\- i9 (— TT < < -k). The superscript I (II) indicates 
the 1st (2nd) sheet [||. 

Here we mention briefly the determination of the un- 
known parameters (d'^-^, Qa and fha) in each model. In 
Models A, B and C, gp is fixed to be 5p/47r = 2.72 for 
reproducing the p-width. In Model B and C, g^ is as- 
sumed to be the same as gp following Ref. |jl^. Re- 
maining parameters d'n and fha are adjusted by mak- 
ing a global fit to the experimental phase shift in the 
I=J=l channel up to 1 GeV Then, d'^^ and d'^Q 

are determined uniquely from d'n by the constraints ob- 
tained from the matching to the 0(p*) chiral lagrangian 
for small s 25 . For Model D, d^j is directly ex- 
tracted from the coefficients L\2{l^) in the O(p^) chiral 
lagrangian. The I—J=l phase shift is not reproduced in 
Model D, since only £2 is considered |2^. The results 
are listed in TABLE || for each model. We take the stan- 
dard choice p=0.77 GeV in determining the parameters. 
We have checked that the experimental phase shift in 
the (/,J)=(0,0) and (/,J)=(2,0) channels is reproduced 
quite well within the experimental error-bars in Model 
B and C below 0.9 GeV. Model A and D also show a 
qualitatively reasonable fit in the above channels at low 
energies. However, Model A underestimates attraction 
for (/,J)=(0,0) above 0.6 GeV. Also, Model D overes- 
timates the attraction and repulsion above 0.4 GeV for 
(/,J)=(0,0) and (/, J)=(2,0), respectively. 

Before making numerical analysis of Djj(s)=0, let us 
first discuss its solution for small values of s where eq. (||) 
may be approximated as 



Djjis) ~ l + (5,,/7r)sln(-s/M2), 



(7) 



in all four models. Here we have neglected 0(s"-^) 
terms and M is a ^-independent scale defined as = 
^2gd (/i)-!r/() Then a solution on the 2nd Riemann sheet 
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FIG. 1; The shift of the pole position in the I—J=l channel 
associated with the decrease of /*. Two kinds of flows exist: 
one toward the origin and the other toward fhp on the real 
axis. Crosses are the pole positions in the vacuum. 



FIG. 2: The shift of the pole position in the I—J=0 channel 
associated with the decrease of /*. Two kinds of flows exist: 
one toward the origin and the other toward fhc on the real 
axis. Crosses are the pole positions in the vacuum. 



in (/,J) = (0,0), (1,1) channels reads 



[W^,{7r/b,,M')\ 

_p*2 



(8) 
(9) 



where W i{z) is the (— iV-th branch of the Lambert W 
function ||2^. From eq.(||) to eq.(^, we have used an 
asymptotic expansion of W-i{z) valid for small and pos- 
itive z !§. Also, F;^ = 4^/; and F* = \/6i^* . Eq.(|) 
explicitly shows that resonance poles appear both in the 
p and the a channels, and are softened in tandem as /* 
decreases in nuclear medium. 

For /* ^ fjr, the small s approximation cannot be jus- 
tified and more complicated behavior of the poles arises. 
To keep track of the trajectory of the poles for wide 
range of /*, we show numerical solutions of eq.(||) for 
the I=J—\ channel (Fig. |l|) and for the I=J=0 chan- 
nel (Fig. ^). In these figures, the crosses indicate the 
position of the poles in the vacuum (/*=/^). A com- 
mon feature in Models A, B and C in the vacuum is that 
there always exists a narrow p {y/s^ — 769 — 76i MeV) 
and a broad/low-mass cr, no matter whether bare reso- 
nance are introduced or not. This is a kind of bootstrap 
situation discussed in [^. Note also that, when the bare 
resonance is introduced, two complex poles appear as a 
result of an interplay between the bare resonance and the 
dynamically generated one. This can be seen in Fig. |l| 
for Models A and C, and in Fig. ^ for Models B and C. 
In Model D where bare resonances are not introduced, 
a broad and low-mass a is dynamically generated while 
narrow p does not appear. 

As /* decreases from its vacuum value, we find two 



types of trajectories in Figs. those moving toward the 
origin and those moving toward rhp(^a-) on the real axis. 
The former trajectories in I—J—1 and J=J=0 channels 
are model independent and correspond to the universal 
softening evaluated in Eq.(^). Since the tt-tt scattering 
amplitude is mainly affected by the poles close to the real 
axis, such soft mode gives a dominant contribution to the 
low-energy amplitude for f*//-^ ^ 1. 

On the other hand, when the system is close to the vac- 
uum state (/t^Z/tt ~ 1), the above universality between p 
and a breaks down. In fact, in the I=J~1 channel, the 
narrow p resonance plays a dominant role and would-be 
soft mode is far away from the real-axis. In the /= J=0 
channel, on the contrary, the broad and low-mass a keeps 
playing a role of the soft mode all the time. 

The similarity and difference between the two channels 
for different values of /* can be seen also by comparing 
the in-medium tt-tt cross section calculated from ajj{s). 

The upper panel of Fig.^ (Fig.^) represents the in- 
medium cross section in Model B for O.S/tt < < J-k 
in the I=J=1 (/= J=0) channel. We find that there is a 
moderate softening -I- broadening of the p-resonance at 
/* = O.T/tt, which is itself an interesting behavior in re- 
lation to the dilepton data in |^. However, even more 
significant change can be seen for the cr-resonance at the 
same value of /*. 

The lower panel of Fig.H (Fig.^) represents the in- 
medium cross section in Model B for 0.1/^ < < 0.5/^ 
in the I=J=1 (/=J=0) channel. We find that sharp 
peaks develop in both channels as /* 0. They are the 
direct consequence of the soft modes located near the ori- 
gin for small /* in Figs. |^||. Namely, both I=J=1 and 
/= J=0 channels are good probes of chiral restoration in 
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FIG. 3: The in-medium tt-tt cross section in the I—J—1 
channel for Model-B. The upper (lower) panel shows the case 
for 0.5/. <f:<fn (0.1/. < /: < 0.5/.). 



the medium if the chiral restoration is substantial. Al- 
though we have taken Model B as an example in Fig. 3-4, 
we have checked that qualitative conclusions are the same 
for other models. 

In summary, we have studied complex poles of the 
in-medium tt-tt scattering amplitude in the I~J~Q and 
I=J=1 channels to explore possible relations between a 
and p in nuclear matter. The N/D method is applied 
to four types of chiral models and the chiral restoration 
is modeled by the decrease of /*. We have found uni- 
versal complex poles in both channels moving in tandem 
toward the origin if /* is sufficiently small. On the other 



FIG. 4: Same with Fig. g for the /= J=0 channel for Model- 
B. 

hand, if /* is not far from its vacuum value, interesting 
non-universal behavior arises and the two resonances act 
rather differently. 

Inclusion of the finite pion mass, the medium effect be- 
yond the mean-field approximation (such as the coupling 
to the particle-hole excitations of the Fermi-sea and the 
difference of /* ^ and / * j, mentioned before) are future 
problems to be examined. Possible connection to other 
theoretical approaches should be also studied. 
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